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ABSTRACT 

We report the discovery of four transiting extrasolar planets (HAT-P-34b-HAT-P-37b) with masses 
ranging from 1.05 to 3.33 Afj and periods from 1.33 to 5.45 days. These planets orbit relatively bright 
F and G dwarf stars (from V = 10.16 to V = 13.2). Of particular interest is HAT-P-34b which is 
moderately massive (3.33 Mj), has a high eccentricity of e = 0.441±0.032 at P = 5.452654±0.000016 d 
period, and shows hints of an outer component. The other three planets have properties that are 
typical of hot Jupiters. 

Subject headings: planetary systems — stars: individual (HAT-P-34, GSC 1622-01261, HAT-P- 
35, GSC 0203-01079, HAT-P-36, GSC 3020-02221, HAT-P-37, GSC 3553-00723 ) 
techniques — spectroscopic, photometric 



1. INTRODUCTION 

Transiting extrasolar planets (TEPs) provide unique 
opportunities to study the properties of planetary ob- 
jects outside of the Solar System. To date well over 100 
such planets have been discovered and characterized 1 "^, 
leading to numerous insights into the physical pro perties 
of pla netary systems (e.g. see the recent review bv lRauerl 
2011). In addition, over a thousa nd strong candidate s 
from Kepler have been identified (jBorucki et al.l 120 111 ), 
greatly expanding our understanding of several aspects of 
planetar y systems, such as the properties of multi-p lanet 
systems ([Latham et a.1.1 [20Tl1 : iLissauer et al]|2011[ ). and 
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the distribution of planetary radii (iRo ward et al.ll2011[) . 
However, due to the large number of important vari- 
ables which influence the physical properties of a planet 
(e.g. its mass, composition, age, irradiation, and tides, to 
name a few), we are still far from an empirically tested, 
comprehensive understanding of the formation and evo- 
lution of planetary systems. 

Here we present the discovery of four new TEPs iden- 
tified by the Hun garian-made Autom ated Telescope Net- 
work (HATNet; IBakos et al.l [200l survey which con- 
tribute to the rapidly-growing sample of TEPs. These 
planets transit relatively bright stars facilitating detailed 
characterization of their properties, such as measure- 
ments of their masses via radial velocity (RV) observa- 
tions of the host star, or measuring their orbital tilt via 
the Rossiter-McLaughlin effect. 

The HATNet survey for TEPs around bright stars 
(9 ^5 r ^5 14.5) operates six wide-field instruments: four 
at the Fred Lawrence Whipple Observatory (FLWO) in 
Arizona (HAT-5, -6, -7, and -10), and two on the roof of 
the hangar servicing the Smithsonian Astrophysical Ob- 
servatory's Submillimcter Array, in Hawaii (HAT-8 and 
-9). Since 200 6, HATNet has anno unced and published 
33 TEPs (e.g. I Johnson et all 120 111 ). In this work we re- 
port our thirty-fourth through thirty-seventh discover- 
ies, around the stars GSC 1622-01261, GSC 0203-01079, 
GSC 3020-02221, and GSC 3553-00723. 

In Section [2] we summarize the detection of the photo- 
metric transit signals and the subsequent spectroscopic 
and photometric observations of each star to confirm the 
planets. In Section |3] we analyze the data to determine 
the stellar and planetary parameters. The properties of 
these planets are briefly discussed in Section H) 

2. OBSERVATIONS 

The observational procedure employed by HATNet 
to discover TEPs has been describe d in detail in sev- 
eral previous discov ery papers (e.g. IBakos et al.l 120101 : 
lLatham et al.l I2009D . In the following subsections we 
highlight specific details of this procedure that are per- 
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TABLE 1 

Summary of discovery data. 



Planet Host GSC 2MASS RA DEC V a Depth b Period 

HH:MM:SS DD:MM:SS mag mmag days 

HAT-P-34 1622-01261 20124688+1806175 20 h 12 m 46.80s +18°06'17.5" 10.162 ±0.073 7.9 5.4527 

HAT-P-35 0203-01079 08130018+0447132 08 h 13 m 00.19s +04°47'13.3" 12.46 + 0.11 9.0 3.6467 

HAT-P-36 3020-02221 12330390+4454552 12 h 33 m 03.96s +44°54'55.3" 12.262 + 0.068 14.7 1.3273 

HAT-P-37 3553-00723 18571105+5116088 18 h 57 m 11.16s +51°16'08.9" 13.23 ± 0.32 c 18.1 2.7974 

a From lDroege et aljfeooa 

Note that the apparent depth of the HATNet transit for all four targets is shallower than the true transit depth due to blending 
with unresolved neighbors in the low spatial resolution HATNet images (the median full-width at half maximum of the point-spread 
function at the center of a HATNet image is ~ 25 ). Also, we applied the trend filtering procedure in non signal-reconstructive mode, 
which reduces the transit depth while increasing the signal-to-noisc ratio of the detection. For each system the ratio of the planet and 
stellar radii, which is related to the true transit depth, is determined in Scction l3.2l using the higher spatial-resolution photometric 
follow-up observations described in Scction l2.4l 
c From lLasker et alj|2008l 



tinent to the discoveries of the four planets presented in 
this paper. 

2.1. Photometric detection 

Table [5] summarizes the HATNet discovery observa- 
tions of each new planetary system. The HATNet images 
were processed and reduced to trend- filtered light curves 
following the pr ocedure described by lBakos et al.1 ()2010[ ) 
and [Pal (20093). The light curves were searched for pe- 
riodic bo x-shaped signals us ing the Box Least-Squares 
(BLS; see lKovacs et al.ll2002D method. We detected sig- 
nificant signals in the light curves of the stars summa- 
rized in Table [1] 

2.2. Reconnaissance Spectroscopy 

High-resolution, low-S/N "reconnaissance" spectra 
were obtained for HAT-P-34 and HAT-P-35 using 
the Tillingh ast Reflector Echelle Spectrograph (TRES; 
lFuresd2008l ) on the 1.5 m Tillinghast Reflector at FLWO. 
These observations were reduc ed and analyzed fo llow- 
ing the procedu r e desc ribed by IQuinn et al.l (|2010l) and 
Buch have et al.l ([2010); the results are listed in Ta- 
ble |3[ For both objects the spectra were single-lined, 
and showed radial velocity (RV) variations on the order 
of ~ 100 ms -1 . Proper phasing of the RV with the pho- 
tometric ephemeris gives confidence in acquiring further, 
high signal-to-noise spectroscopic observations to refine 
the orbit (see §l2~3"]). While for HAT-P-34 the variations 
initially did not appear to phase with the photometric 
ephemeris, we entertained the possibility of a very signif- 
icant non-zero eccentricity, and pursued follow-up of the 
target. For HAT-P-35 the variations were in phase with 
the photometric ephemeris indicating a ~ 2.7 Mj com- 
panion. For both HAT-P-36 and HAT-P-37 we obtained 
two TRES spectra near each of the predicted quadra- 
ture phases. For both objects the spectra were single- 
lined. For HAT-P-36 the resulting RV measurements 
showed ~ 400 ms -1 variation in phase with the pho- 
tometric ephemeris, while for HAT-P-37 the RV mea- 
surements showed ~ 260 ms -1 variation in phase with 
the ephemeris. We opted to continue observing both of 
these objects using TRES with the aim of confirming the 
planets. The TRES observations of HAT-P-36 and HAT- 
P-37 are discussed further in the following subsection. 

2.3. High resolution, high S/N spectroscopy 



We proceeded with the follow-up of each candidate by 
obtaining high- resolution, high-S /N spectra to character- 
ize the RV variations, and to refine the determination of 
the stellar parameters. These observations are summa- 
rized in Table |H The RV measurements and uncertain- 
ties for HAT-P-34 through HAT-P-37 are given in the 
Appendix, and in Tables [71-HTJ1 respectively. The period- 
folded data, along with our best fit described below in 
Section [3l are displayed in Figures [SHU 

Four facilities were used in the confirmation of these 
planets (including three separate fa cilities used for H AT- 
P-34). These facilities are HIRES (jVogt et al J 1 1991 on 
the 10 m Keck I teles cope in Hawaii, the H igh-Dispersion 
Spectrograph (HDS; iNoguchi et all 120021 ) on the 8.3 m 

Subaru telescope in Hawaii, the Fibre-fed Echelle Spec- 
trograph (FIES) on the 2.5 m Nordic Optical Telescop e 
(NOT) at La Palma, Spain (|Diupvik fc Anderseril2lH(3 l. 
and TRES on the FLWO 1.5 m telescope. 

The HIRES and HDS observations made use of the 
iodine-cell method (jMarcv fc Butlerl 119921: IButler et all 
1996) for precise wavelength calibration and relative 
RV determination, while the FIES and TRES obser- 
vations made use of Th-Ar lamp spectra obtained be- 
fore and after the science exposures. The HIRES ob- 
servations were redu ced to relativ e RVs i n the barycen- 
tric fr ame f ollowing IButler et all (|1996l) . Uohnson et al.l 
(|2009fl . and lHoward et al.1 (120111. t he HPS ob servations 
were reduced following ISato et all (pfM 120051 ). and the 
FIES and TRES obs ervations were reduced following 
iBuchhave et~aH (|2010t ). 

We found that for all four systems the RV residuals 
from the best-fit models, described below in Section l3~2l 
exhibit excess scatter over what is expected based on the 
formal measurement uncertainties. Such excess scatter, 
or "jitter" has been well known for stars, and can stem 
for multiple sources. The excess is in the residuals of 
the observations with respect to a physical (and possibly 
instrumental) model. If this model is not adequate, the 
residuals can be larger than expected. For example, in 
the case of HAT-P-34b, ignoring the linear trend in the 
RVs would lead to a much increased "jitter" . Additional 
planets may cause jitter, as the limited number of RV ob- 
servations is not enough to uniquely identify and model 
such systems. The typical source of the jitter, however, 
is the star itself, namely inhomogeneitie s (spots, flares, 
plages, etc.) on the stellar surfac e (e.g. iMakarov et al.l 
120091 : iMartmez-Arnaiz et al.l 1201(3 ) causing jitters up to 
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TABLE 2 

Summary of photometric observations 



Instrument /Field 


Date(s) 


Number of Images 


Cadence (sec) 


Filter 


HAT-P-34 










HAT-7/G293 


2008 Oct-2009 May 


755 


330 


T 


HAT-8/G293 


2008 Sep-2008 Dec 


2611 


330 


r 


HAT-6/G341 


2007 Sep-2007 Dec 


1949 


330 


R 


HAT-9/G341 


2007 Sep-2007 Nov 


2379 


330 


R 


KeplerCam 


2010 May 21 


263 


60 


z 


KeplerCam 


2010 Oct 10 


530 


30 


i 


II A 1 I > qc 
1 1 ; l I -r-oo 










HAT-5/G364 


2009 May 


21 


330 


V 


HAT-9/G364 


2008 Dec-2009 May 


3155 


330 


r 


KeplerCam 


2011 Jan 16 


110 


100 


i 


FTN 


2011 Jan 23 


185 


45 


i 


KeplerCam 


2011 Mar 08 


268 


60 


i 


HAT-P-36 










HAT-5/G143 


2010 Apr-2010 Jul 


4471 


210 


r 


HAT-8/G143 


2010 Apr-2010 Jul 


6262 


210 


r 


KeplerCam 


2010 Dec 24 


131 


130 


i 


KeplerCam 


2011 Feb 03 


101 


100 


i 


KeplerCam 


2011 Feb 07 


105 


100 


i 


KeplerCam 


2011 Feb 15 


186 


60 


i 


HAT-P-37 










HAT-7/G115 


2009 Sep-2010 Jul 


7102 


210 


r 


HAT-9/G115 


2008 Aug-2008 Sep 


2293 


330 


R 


KeplerCam 


2011 Feb 23 


37 


165 


i 


KeplerCam 


2011 Mar 23 


73 


134 


i 


KeplerCam 


2011 Apr 06 


102 


134 


i 



TABLE 3 

Summary of reconnaissance spectroscopy 
observations 3, 



Instrument 


HJD 


7RV b 
(kms -1 ) 


CC Peak c 


HAT-P-34 








TRES 


2454935.00839 


-47.81 


0.726 


TRES 


2454966.97204 


-49.19 


0.939 


TRES 


2454998.97956 


-49.12 


0.819 


HAT-P-35 








TRES 


2455289.64284 


41.24 


0.863 


TRES 


2455291.62482 


40.54 


0.883 


TRES 


2455320.64182 


40.83 


0.935 


TRES 


2455321.64383 


41.08 


0.940 



a For HAT-P-36 and HAT-P-37, which were confirmed us- 
ing the TRES spectrograph, there is no clear distinction be- 
tween reconnaissance and high-precision observations. We 
do not list the results from the analysis of the TRES spec- 
tra for these targets here, these are instead described in 
Section nOl 

The heliocentric RV of the target in the IAU system, 
and corrected for the orbital motion of the planet. 
c The peak value of the cross-correlation function be- 
tween the observed spectrum and the best-matching syn- 
thetic template spectrum (normalized to be between and 
1). Observations with a peak height closer to 1.0 generally 
correspond to higher S/N spectra. 
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Fig. 1. — HATNet light curves of HAT-P-34 through HAT-P-37. See Table [2] for a summary of the observations. For each planet we 
show two panels. The top panel shows the unbinned light cu rve f olded with the period resulting from the global fit described in Section [3] 
The solid line shows the model fit to the light curve (Section l3.2l l. The bottom panel shows the region zoomed-in on the transit. The dark 
filled circles show the light curve binned in phase with a binsize of 0.002. 
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100 ms. Granulation and stellar oscillations contribute 
at a smaller scale, but are present for non- active stars 
that ar e outside the insta bility strip. A recent publica- 
tion bv lCegla et al.1 (j2012f ) discusses the stellar jitter due 
to variable gravitational redshift of the star, as the stellar 
radius changes due to oscillations (AR of 10 -4 causing 
~ 0.1ms" 1 ). And, of course, systematics in the instru- 
ment further inflate the jitter. A review of RV jitter of 
stars o bserved by the Keck telescope are given in Wright 
(l200l. 

In order to ensure realistic estimates of the system pa- 
rameter uncertainties we add in quadrature an RV jitter 
to the formal RV measurement uncertainties such that x 2 
per degree of freedom is unity for the best-fit model for 
each planet. We adopt an independent jitter for the ob- 
servations made by each instrument of each planet. The 
RV uncertainties given in Tables [THTOl do not include this 
jitter; we do include the jitter in Figures [2H5] 

For HAT-P-34 and HAT-P-35 we also show the S in- 
dex, which is a measure of the chromospheric activity 
of the star derived from the flux in the cores of the 
Ca II H and K lines. This in dex was computed follow- 
ing llsaacson fc Fischer! (12010ft and has been calibr ated 
to the scale of iVaughan, Preston fc Wilson! (|1978f ). A 
procedure for obtaining calibrated S index values from 
the TRES spectra has not yet been developed, so we 
do not provide these measurements for HAT-P-36 or 
HAT-P-3 7. We convert the S index values to logi?^ K 
following iNoves et al.l (|1984l ) and find median values of 
logi^K = -4.859 and logi^K = -5.242 for HAT-P- 
34 and HAT-P-35, respectively. These values imply that 
neither star has a particularly high level of chromospheric 
activity. 

Following lQueloz et al.l (|200lD and lTorres et al.l (|2007|) . 
we checked whether the measured radial velocities are 
not real, but are instead caused by distortions in the 
spectral line profiles due to contamination from a nearby 
unresolved eclipsing binary. A bisector (BS) analysis 
for each system based on the K eck and TRES spec - 
tra was done as described in §5 of iBakos et al.l (|2007aD . 
For HAT-P-35, which is relatively faint, we found that 
the measured BSs were significantly affected by scat- 
tered moonlight and applied an emp i rical correction 
for this effect follow ing lHartman et al.1 (|2009f ) (see also 
iKovacs et all (|2010f )). For HAT-P-34 the BS scatter is 
fairly high (~ 25 ms -1 ), but this is in line with the 
high RV jitter (~ 60ms -1 ), which is typical of an F 
star with ^sinj = 24 .0 ± 0.5 km s" 1 (|Saar et al.l 120031 : 
lHartman et~al]|2011b| ). 

None of the systems show significant bisector span vari- 
ations (relative to the semi- amplitude of the RV varia- 
tions) that phase with the photometric ephemeris. Such 
variations are generally expected if the transit and RV 
signals were due to blends rather than planets. While the 
lack of bisector span variations does not exclude all blend 
scenarios, it does significantly limit the possible blend 
scenarios that can reproduce our current data within the 
measurement errors, i.e. configurations that are compati- 
ble with the photometric and spectroscopic observations, 
proper motions, color indices, and moderately high res- 
olution imaging. We have found in the past that invok- 
ing detailed blend modeling to exclude all possible blend 
configurations and confirm the planet hypothesis (e.g. 
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Fig. 2. — Top panel: High-precision RV measurements for 
HAT-P-34 shown as a function of orbital phase, along with 
our best-fit model (see Table [BJ. Open triangles show mea- 
surements from Subaru/HDS, filled circles show measurements 
from Keck/HIRES, and filled triangles show measurements from 
NOT/FIES. Zero phase corresponds to the time of mid-transit. 
The center-of-mass velocity and a linear trend have been sub- 
tracted. Second panel: velocity O — C residuals from the best-fit 
single Keplerian orbit model as a function of time. The residuals 
show a slight linear trend, possibly indicating a third body in the 
system. Note that the zero-points of the three separate instru- 
ments are independently free parameters. Third panel: velocity 
O — C residuals from the best fit including both the Keplerian or- 
bit and linear trend, shown as a function of orbital phase. The 
error bars include a jitter term (56.0ms _1 for the Keck/HIRES 
observations, and 32.0ms _ for the Subaru/HDS observations; no 
jitter has been added to the NOT/FIES RV u ncer tainties) added 
in quadrature to the formal errors (see Section 13.211 . Fourth panel: 
bisector spans (BS) from Kcck/HlRES, with the mean value sub- 
tracted. The measurement from the template spectrum is included. 
Bottom panel: Chromospheric activity index S measured from the 
Keck spectra. Note the different vertical scales of the panels. Ob- 
servations shown twice are represented with open symbols. 

lHartman et al.ll2011al sections 3.2.2, 3.2.3) is rarely of 
any incremental value when the ingress and egress dura- 
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TABLE 4 

Summary of high-SN spectroscopic observations 
used in measuring the orbits 



Instrument 



Date(s) 



Number of 
RV obs. 



HAT-P-34 

Subaru/HDS 
Keck/HIRES 
NOT/FIES 

HAT-P-35 

Keck/HIRES 
NOT/FIES 

HAT-P-36 



2010 May 

2010 Jun-2010 Sep 

2010 Jul-2010 Aug 



2010 Scp-2010 Dec 
2010 Oct 



6 

14 
10 



FLWO 1.5/TRES 2010 Dec-2011 Jan 



HAT-P-37 

FLWO 1.5/TRES 



2011 Mar-2011 May 



12 



13 



a One of the NOT/FIES spectra of HAT-P-35 was aborted 
early due to morning twilight and high humidity, another ex- 
posure was obtained partly during transit and may be af- 
fected by the Rossitcr-McLaughlin effect. The remaining three 
NOT/FIES spectra do not provide sufficient phase coverage to 
constrain the orbit. We therefore do not include the velocities 
measured from these spectra in the analysis of HAT-P-35. 
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Fig. 3.— Keck/HIRES observations of HAT-P-35. The panels 
are as in Figure [2] The parameters used in the best-fit model are 
given in Table [6] 



tions are short relative to the total transit duration, the 
RV variations exhibit a Keplerian orbit in phase with the 
photometric ephemeris, and bisector spans show no cor- 
relation with the orbit. We conclude that the velocity 
variations detected for all four stars are real, and that 
each star is orbited by a close-in giant planet. 




0.4 0.6 0.8 

Phase with respect to T c 

Fig. 4.— FLWO 1.5m/TRES observations of HAT-P-36. The 
panels are as in Figure [2] The S index is not available for these 
observations. The parameters used in the best-fit model are given 
in Table El 

2.4. Photometric follow-up observations 

In order to permit a more accurate modeling of the 
light curves, we conducted additional photometric ob- 
servations using the KeplerCam CCD camera on the 
FLWO 1.2 m telescope, and the Spectral Instrument 
CCD on the 2.0 m Faulkes Telescope North (FTN) 
at Haleakala Observatory in Hawaii, which is oper- 
ated by the Las Cumbres Observatory Global Telescope 
(LCOGT). The observations for each target are summa- 
rized in Table H 

The reduction of these images was performed as de- 



HAT-P-34b, HAT-P-35b, HAT-P-36b, and HAT-P-37b 



7 




Phase with respect to T c 

Fig. 5.— FLWO 1.5 m/TRES observations of HAT-P-37. The 
panels are as in Figure [2] The S index is not available for these 
observations. The parameters used in the best-fit model are given 
in Table [6] 
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Fig. 6. — Unbinned transit light curves for HAT-P-34, acquired 
with KeplerCam at the FLWO 1.2 m telescope. The li ght c urves 
have been EPD- and TFA-processed, as described in § 13.21 The 
dates of the events are indicated. Curves after the first are dis- 
placed vertically for clar ity. Our best fit from the global modeling 
described in Section l3.2l is shown by the solid lines. Residuals from 
the fits are displayed at the bottom, in the same order as the top 
curves. The error bars represent the photon and background shot 
noise, plus the readout noise. 
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Fig. 7. — Similar to Figure [6] here we show the follow-up light 
curves for HAT-P-35. 
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Fig. 8. — Similar to Figure(6] here we show the follow-up light 
curves for HAT-P-36. 



scribed by iBakos et al.l (|2010D . We applied Exte rnal Pa- 
rameter Decorrelation (EPD; lBakos et aLll2010f) and the 
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Fig. 9. — Similar to Figure [6] here we show the follow-up light 
curves for HAT-P-37. 



Trend Filtering Algorithm fTFA: lKovacs et al]l200l to 
remove trends simultaneously with the light curve mod- 
eling. The final time series, together with our best-fit 
transit light curve models, are shown in the top por- 
tion of Figures MM for HAT-P-34 through HAT-P-37, 
respectively. The individual measurements, permitting 
independent analysis by other researchers, are reported 
in the Appendix, in Tables [iTrfMl (the full data are avail- 
able in electronic format). 

3. ANALYSIS 
3.1. Properties of the parent star 

Stellar atmospheric parameters for HAT-P-34 and 
HAT-P-35 were measured using our template spectra ob- 
tained with the Keck/HIRES instrument, and the analy- 
sis package known as Sp ectroscopy Made Easy (SME; 
iValenti &: Piskunovl Il996fl. alo n g wit h the atomic line 
database of iValenti fc Fischer! (]2005h . For HAT-P-36 
and HAT-P-37 the stellar atmospheric parameters were 
determined by cross-correlating the TRES observations 
ag ainst a finely s ampled grid of synthetic spectra based 
on iKuruczl (|2005l ) model atmospheres. This procedure, 
known as Stellar Parameter Classification (SPC), will be 
described in detail in a forthcoming paper (Buchhave 
et al., in preparation). We note that SPC has been 
performed in the past on numerous HATNet transiting 
planet candidates (Buchhave, personal communication), 
and the results were consistent with those of SME. 

For each star, we obtained the following initial spec- 
troscopic parameters and uncertainties: 

• HAT-P-34 ~ effective temperature T cS + = 6400 ± 
100 K, metallicity [Fe/H] = +0.21 ± 0.1 dex, stellar 



surface gravity logg* = 3.98 ± 0.1 (cgs), and pro- 
jected rotational velocity v sin i = 24.5±1.0kms~ 1 . 

• HAT-P-35 - effective temperature T cff * = 5940 ± 
88 K, metallicity [Fe/H] = +0.01 ± 0.08 dex, stellar 
surface gravity logg* = 3.98 ± 0.1 (cgs), and pro- 
jected rotational velocity vsmi = 0.5 ±0.5 km s -1 . 

• HAT-P-36 - effective temperature T oSir = 5850 ± 
100 K, metallicity [Fe/H] = +0.38 ± 0.1 dex, stellar 
surface gravity logg* = 4.73 ± 0.17 (cgs), and pro- 
jected rotational velocity v sin i = 2.86+0.5 km s _1 . 

• HAT-P-37 - effective temperature T cS * = 5570 ± 
100 K, metallicity [Fe/H] = +0.09 ± 0.1 dex, stellar 
surface gravity logg* = 4.67 ± 0.1 (cgs), and pro- 
jected rotational velocity v sin i = 2.95+0.5 km s . 



Following iBakos et al.l ((2010), these initial values 
of T^ff*, log<?*, and [Fe/H] were used to determine 
the quadratic limb-darkening coefficients needed in the 
global modeling of the follow-up photometry (summa- 
rized in Section l3.2p . This analysis yields p*, the mean 
stellar density, which is closely related to a/R+, the nor- 
malized semimajor axis, and provides a tighter constraint 
on the stellar parameter s than does the spect roscopically 
determined logg* (e.g. S ozzetti et al.l l2007h . We com- 
bined /?+, Toff*, and [Fe/H] with stell ar evoluti o n mod els 
from the Yonsei-Yale (YY) series by lYi et al.l (|2001h to 
determine probability distributions of other stellar prop- 
erties, including log<?*. For each system we carried out 
a second SME or SPC iteration in which we adopted 
the new value of logg* so determined and held it fixed 
in a new SME or SPC analysis, adjusting only T e g+, 
[Fe/H], and vsmi, followed by a second global modeling 
of the RV and light curves, together with improved limb 
darkening parameters. The final atmospheric parameters 
that we adopt, together with stellar parameters inferred 
from the YY models (such as the mass, radius and age) 
are listed in Table [5] for all four stars. 

The inferred location of each star in a diagram of a/i?* 
versus Tcff*, analogous to the classical H-R diagram, is 
shown in Figure 1101 In each case the stellar proper- 
ties and their la and 2a confidence ellipses are displayed 
against the backdrop of model isochrones for a range of 
ages, and the appropriate stellar metallicity. For compar- 
ison, the locations implied by the initial SME and SPC 
results are also shown (in each case with a triangle). 

The stellar evolution modeling provides color indices 
that we compared against the measured values, as a san- 
ity check. For each star we used t he near-infrared magn i- 
tudes from the 2MASS Catalogue (jSkrutskie et alJ200fl ). 
which are given in Table [5l These were converted to the 
photometric s ystem of the mode ls (ESO) using the trans- 
formations bv lCarpenteH (|200lD . The resulting 2MASS- 
based color indices were all consistent (within la) with 
the stellar model based color indices. 

The distance for each star given in Table [5] was com- 
puted from the absolute K magnitude from the models 
and the 2MASS K s magnitudes, ignoring extinction. 

3.2. Global modeling of the data 

We modeled simultaneously the HATNet photome- 
try, the follow-up photometry, and the high-precision 
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Fig. 10. — Comparison of the measured values of T e ff* a nd a/FU for HA T-P-34 (upper left), HAT-P-35 (upper right), HAT-P-36 (lower 
left) and HAT-P-37 (lower right) to model isochrones from Yi ct al. (2001). The isochrones are generated for the measured metallicity of 
each star, and for ages of 0.5 Gyr and 1 to 3 Gyr in steps of 0.25 Gyr for HAT-P-34, and of 0.5 Gyr and 1 to 14 Gyr in steps of 1 Gyr for 
HAT-P-35, HAT-P-36 and HAT-P-37 (ages increase from left to right in each plot). The lines show the la and 2<r confidence ellipses for 
the measured parameters. The initial values of T e ff+ and a/ij* from the initial spectroscopic and light curve analyses are represented with 
a triangle in each panel. 



RV measurements using the procedures described by 
Bak os et all (|2010f ). Namely, the best fit was determined 
by a downhill simplex minimization, and was followed 
by a Monte-Carlo Markov Chain run to scan the param- 
eter space a round the minimum, and est ablish the errors 
iPlll (I2009bl) . For each system we used a lMandel & Agoll 
( 2002) transit model, together with the EPD and TFA 
trend- filters, to descri be the follow-up light curves, a 
IMandel fc Agoll (|2002D transit model for the HATNet 
lig ht curve (s), and a Keplerian orbit using the formalism 
of EH (|2Q09ah for the RV curve(s). For HAT-P-34 we 
included a linear trend in the RV model, but find that 
it is only significant at the ~ 2a level; the planet and 
stellar parameters are changed by less than la when the 
trend is not included in the fit. The parameters that we 
adopt for each system are listed in Table |U In all cases 
we allow the eccentricity to vary so that the uncertainty 
on this parameter is propagated into the uncertainties 
on the other physical parameters, such as the stellar and 
planetary masses and radii; the observations of HAT- 
P-35b, HAT-P-36b, and HAT-P-37b are consistent with 
these planets being on circular orbits. 

4. DISCUSSION 

We have presented the discovery of four new transiting 
planets. Below we briefly discuss their properties. 

4.1. HAT-P-34b 



HAT-P-34b is a relatively massive M p = 3.328 ± 
0.211 Mj planet on a relatively long period (P = 
5.452654 ± 0.000016 d), eccentric (e = 0.441 ± 0.032) or- 
bit. There are only five known transiting planets with 
higher eccentric i ties (HAT-P-2b, e = 0.5171 ± 0.0033, 
iPal et alJ [2010t iBakos et alj I2007at CoRoT-lOb, e = 
0.53± q.04.lBonomo et al.ll2010t CoRoT-20b. e = 0.562± 
0.013, iDeleuil et alJ 120121: HP 171 56b, e = 0.669 ± 
0.008, IMadhusudhan fc Winnl 12009) and HD 80606b, 
e = 0.9330 ± 0.0005. IHe^rard'etlifllMch . all of which 
have longer orbital periods than HAT-P-34b. Of these 
planets, HAT-P-2b is most similar in orbital period to 
HAT-P-34b, but it has a mass that is more than two 
times larger than that of HAT-P-34b. Two planets with 
masses, radii and equilibrium temperatures within 10% 
of the values of HAT-P-34b (assuming z ero albedo and 



full h eat redistributio n) are CoRoT-18b Hcbrar d et al.1 
1201 If ) and WASP-32b (jMaxted et al.ll2010f) : however nei- 
ther of these planets has a significant eccentricity. 

HAT-P-34b is a promisin g target for measuring the 
Rossiter-McLaughlin effect (|Rossiterl Il924t iMcLaughlir] 
11924f t. since the host star is bright (V = 10.16), has 
a significant spin (usinz= 24.0 ± O.Skms^ 1 ), and the 
transit is moderately long (T14 = 0.1455 ± 0.0016 days). 
Also, the transit is far from equatorial (b = 0.336toi2s)> 
a configuration that is important for resolving the degen- 
eracy between vsini and A, which is the sky-plane pro- 
jected angle betwee n the planetary or bital normal and 
the stellar spin axis. I Winn et"al] (|20100 pointed out that 
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TABLE 5 

Stellar parameters for HAT-P-34 through HAT-P-37 





HAT-P-34 


HAT-P-35 


HAT-P-36 


HAT-P-37 




Parameter 


Value 


Value 


Value 


Value 


Source 



Spectroscopic properties 

T cff * (K) 6442 ± 88 

[Fe/H] +0.22 ± 0.04 

tjsini (bus" 1 ) . . . 24.0 ± 0.5 

fmac (kms -1 ) .... 5.05 

u m i c (kms -1 ) .... 0.85 

7rv (kms" 1 ) -49.26±0.30 

Photometric properties 

V (mag) 10.162 ±0.073 

V-I c (mag) 0.557 ±0.12 

J (mag) 9.460 ± 0.022 

H (mag) 9.322 ± 0.030 

K s (mag) 9.247 ± 0.023 

Derived properties 

M+ (M Q ) 1.392 ±0.047 

R* (Re) i-535lg;Jg 

logp* (cgs) 4.21 ±0.06 

L*(Lq) 3.63lg : ™ 

M v (mag) 3.32 ±0.19 

M K (mag.ESO) . . 2.24 ± 0.17 

Age (Gyr) 1.7+°;* 

Distance (pc) 2571^ 



6096 ± 88 
+0.11 ±0.08 
0.5 ±0.5 
4.52 
0.85 
40.95 ± 0.20 



12.46 ±0.11 
0.662 ± 0.12 
11.358 ±0.024 
11.072 ±0.023 
11.030 ±0.021 



1.236 ± 0.048 
1.435 ± 0.084 
4.21 ±0.04 

9 ctr+0.40 
z - oo -0.30 

3.77 + 0.15 
2.43 ±0.13 

o c+0.8 
J - o -0.5 

535 ± 32 



5560 ± 100 
+0.26 ± 0.10 
3.58 ±0.5 
0.00 
0.00 
-16.29 + 0.10 



12.262 ± 0.068 
0.760 ±0.13 
11.046 ±0.027 
10.723 ± 0.030 
10.603 ± 0.021 



1.022 ± 0.049 
1.096 ± 0.056 
4.37 ±0.04 
1.03 ±0.15 
4.83 ±0.17 
3.14 ±0.12 
6.6±?;1 
317 ±17 



5500 ± 100 
+0.03 + 0.10 
3.07 + 0.5 



13.23 + 0.32 

12.092 ± 0.027 
11.714 + 0.032 
11.667 ±0.020 



0.929 ± 0.043 

n O77+0.059 

4.52 + 0.05 
62+ - 11 
5.41 + 0.19 
3.64 + 0.14 

411 + 26 



Spec. Analysis.' 
Spec. Analysis. 
Spec. Analysis. 
Spec. Analysis. 
Spec. Analysis. 



-20.53 + 0.1 TRES 



TASS,GSC b 

TASS 

2MASS 

2MASS 

2MASS 



YY+a/i?*+Spec. 
YY+a/i?*+Spec. 
YY+a/ii^+Spec. 
YY+a/R^+Spec. 
YY+a/i?*+Spec. 
YY+a/iJ*+Spec. 
YY+a/i?*+Spec. 
YY+a/ii^+Spec. 



Analysis. 
Analysis 
Analysis 
Analysis 
Analysis 
Analysis 
Analysis 
Analysis 



a Based on the analysis of high resolution spectra. For HAT-P-34 and HAT-P-35 this corresponds to SME applied to iodine-free 
Keck/HIRES spectra, while for HAT-P-36 and HAT-P-37 this corresponds to SPC applied to the TRES spectra f Section |3~H . These 
parameters also have a small dependence on the iterative analysis incorporating the isochronc search and global modeling of the data, 
as described in the text. 

b For HAT-P-34 through HAT-P-36 the value is taken from the TASS catalog, while for HAT-P-37 the value is taken from the GSC 
version 2.3.2. 

c YY+a/7?* -t-Spcc. Analysis — Based on the YY isochroncs (Yi ct al. 2001), a/R* as a luminosity indicator, and the spectroscopic 
analysis results. 



hot Jupiters around stars with T^* > 6250 K have a 
higher chance of being misaligned. Based on the effec- 
tive temperature of the host star 6442 ± 88 K, we thus 
expect that HAT-P-34b has a higher chance of misalign- 
ment (note that this may not necessarily yield a non- 
zero A, if HAT-P-34b 's orbit is ti l ted al ong the line of 
sight). Alternatively, [Schlaufman (2010) used a stellar 
rotation model and observed v sin i values to statistically 
identify TEP systems that may be misaligned along the 
line of sight, and concluded these preferentially occur at 
M* > 1.2Mq. Based on the stellar mass alone (1.39 M Q ) 
we the chances for misalignment are increased. 

4.2. HAT-P-35b 

HAT-P-35b is a very typical M p = 1.054 ± 0.033 Mj, 
R p = 1.332 ± 0.098 i?j planet on a P = 3.646706 ± 
0.000021 d orbit and with an equilibrium temperature 
of T cq = 1581 ± 45 K (again, assuming zero albedo 
and full heat redistribution). There are four other 
planets with masses, radii and equilibrium tempera- 
tures that are all within 1 0% of the values fo r HAT-P- 
35b. These are HAT-P -5b (|Bakos et al. l l2007bD. HAT-P- 
6b (INoves et al.l 120081) OGLE-TR-211b"(|UdaLski et all 
I2008D, and WASP-26b (ISmallev et al.ll2010D . The stellar 
effective temperature (6096 ± 88 K) is close to the as- 
sumed border-line between well-aligned and misaligned 
systems, making it an interesting system for testing the 
RM effect (with the caveat that wsini, and thus the ex- 
pected amplitude of the anomaly, is low). 



4.3. HAT-P-36b 

HAT-P-36b is a very short period (P = 1.327347 ± 
0.000003d) planet with a mass of M p = 1.832±0.099 Mj, 
a radius of R p — 1.264 ± 0.071 and an equilibrium 
temperature of T eq = 1823 ± 55 K. There are two other 
planets with masses, radii, and equilibrium tempera- 
tures within 10% of the values for HAT-P -36b: TrES-3b 
(lO'Donovan et al.l I2007T ) and WASP-3b (jPollacco et al.1 
2008). 

4.4. HAT-P-37b 

Like the preceding planets, HAT-P-37b also has very 
typical physical properties, with M p = 1.169 ± 0.103 Mj, 
k p = 1.178 ± 0.077 Rj, P = 2.797436 ± 0.000007 d, 
and T oq = 1271 ± 47 K. Three planets with masses, 
radii and equilibrium temperatures withi n 10% of the 
value s for HAT-P-37b a re HP 189733b dBouchv et al.1 
l2005 | ). OGLE-TR-11 3b (jBouchv et aLll2004D . and XO- 
5b (|Burke et al.ll2008h . HAT-P-37 lies just outside of the 
field of view of the Kepler Space mission and is listed in 
the Kepler Input Catalog (KICE3) as KIC 12396036. 

4.5. On the eccentricity of HAT-PS^b 

According to Ad ams and Laughlinl (|2006[ ). the eccen- 
tricity of a hot Jupiter's orbit decays both due to the 
tides on the star and due to the tides on the planet, 

14 http:/ /www. cfa.harvard.edu/kepler/kic/kicindex. html 
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TABLE 6 

Orbital and planetary parameters for HAT-P-34b through HAT-P-37b 





HAT-P-34b 


HAT-P-35b 


HAT-P-36b 


HAT-P-37b 


Parameter 


Value 


Value 


Value 


Value 



Light curve parameters 

P (days) 5.452654 ± 0.000016 

T c (BJD) a 2455431.59629 ± 0.00055 

T w (days) a 0.1455 ±0.0016 

T12 = T 34 (days) a 0.0121 ± 0.0013 

a/R+ 9.48 ±0.64 

C/R* 14.99 ±0.09 

Ftp/ Ft* 0.0801 ± 0.0026 



3.646706 ± 0.000021 
2455578.66081 ± 0.00050 
0.1640 ±0.0018 
0.0162 ±0.0017 
7.45 ±0.37 
13.52 ±0.09 
0.0954 ± 0.0027 



b = a cos i I i?* 

i (deg) 

Quadratic limb-darkening coefficients b 
ci, i (linear term) 



0.113 
0.336 
87.1 ± 1.2 



+0.080 
0.062 
+0.099 
0.128 



0.128 

0.357 
87.3 ± 1.0 



+0.078 
0.066 
+0.092 
0.127 



C2, 
Cl, 
C2, 



(quadratic term) 



RV parameters 

K (ms" 1 ) 

e cos u) c 

e sin u) c 

e 

w (deg) 

7 (ms — 1 d _1 ) 

RV jitter 

Keck/HIRES (ms" 1 ) 

Subaru/HDS (ms" 1 ) 

NOT/FIES (ms" 1 ) 

FLWO 1.5/TRES (ms" 1 ) . 

Secondary eclipse parameters 

T a (BJD) 

T Sj l4 

T S> 12 

Planetary parameters 

M p (Mj) 

Rp 0Rj) 

C{M v ,Rp) d 

p v (gem" 3 ) 

log9p (cgs) 

a (AU) 

T eq (K) 

<->' 

(F) (10 9 ergs- 1 cm- 2 ) f ... 



0.1785 
0.3825 
0.1269 
0.3728 



343.1 ±21.3 
0.410 ±0.031 
0.156 ±0.052 
0.441 ± 0.032 
20 ± 14 
0.8683 ± 0.4719 



56.0 
32.0 
0.0 



2455435.721 ±0.099 
0.1871 ±0.0170 
0.0176 ±0.0052 



3.328 ±0.211 
1 197+0.128 

1 - ly '-0.092 

0.23 
2.40 ± 0.63 
3.76 ± 0.08 
0.0677 ± 0.0008 
1520 ± 60 
0.269 ±0.029 

i 91 +0.23 
1 - z± -0.16 



0.2198 
0.3587 



120.7 ±2.2 
-0.004 ±0.013 
-0.017 ±0.026 

0.025 ±0.018 
248 ± 93 



3.7 



2455580.476 ± 0.030 
0.1596 ±0.0076 
0.0156 ±0.0019 



1.054 ±0.033 
1.332 ±0.098 
0.49 
0.55 ±0.11 
3.17 ±0.06 
0.0498 ± 0.0006 
1581 ±45 
0.064 ± 0.005 

i -^ i -0.14 



1.327347 ± 0.000003 
2455565.18144 ± 0.00020 
0.0923 ± 0.0007 
0.0107 ±0.0007 
4.66 ± 0.22 
24.51 ± 0.14 
0.1186 ± 0.0012 
n nq 7 +0.057 
u.uyi-_ 048 

n oi 9 +0.078 
u.orz_ 105 

86.0 ± 1.3 



0.3142 
0.3113 



334.7 ±14.5 
-0.002 ± 0.032 
0.051 ± 0.040 
0.063 ± 0.032 
95 ±63 



2.797436 ± 0.000007 
2455642.14318 ± 0.00029 
0.0971 ±0.0015 
0.0153 ±0.0013 
9 32+ ' 42 

t, - J -0.57 

24.33 ±0.18 
0.1378 ±0.0030 

n okc+0.044 

5n c i+ ' 041 
"• JUJ -0 .062 



86.9 



+0.4 
0.5 



0.3156 
0.3032 
0.2477 
0.3082 



177.7 ± 14.8 
-0.017 ±0.039 
0.007 ± 0.060 
0.058 ± 0.038 
164 ± 84 



33.6 



2455565.844 ± 0.027 
0.1013 ± 0.0071 
0.0120 ± 0.0015 



1.832 ± 0.099 
1.264 ±0.071 
0.11 
1.12 ± 0.19 
3.45 ± 0.05 
0.0238 ± 0.0004 
1823 ± 55 
0.067 ± 0.005 
2.49 ± 0.30 



25.8 



2455643.512 ± 0.070 
0.0981 ± 0.0083 
0.0153 ±0.0029 



1.169 ±0.103 
1.178 ±0.077 
0.02 
0.89 ±0.19 
3.32 ± 0.07 
0.0379 ± 0.0006 
1271 ± 47 
0.081 ±0.009 

H0. 102 



0.589" 



0.075 



a T c : Reference epoch of mid transit that minimizes the correlation with the orbital period. T14: total transit duration, time between first to last 
contact; T12 = T34: ingress/egress time, time between first and second, or third and fourth contact. Barycentric Julian dates (BJD) throughout 
the paper are calculated from Coordinated Universal Time (UTC). 

b Values for a quadratic law, adopted from the tabulations by Claret (2004) according to the spectroscopic (SME) parameters listed in Table [5] 
c Lagrangian orbital parameters derived from the global modeling, and primarily determined by the RV data. 
d Correlation coefficient between the planetary mass M p and radius Ftp. 

e The Safronov number is given by G = i(V CS c/V orb ) 2 = (a/Rp)(M p /M*) Csee lHansen He Barman|[2007D . 
1 Incoming flux per unit surface area, averaged over the orbit. 



with the tides on the planet dominating the circular- 
ization as long as the tidal quality factor of the planet 
(Qp) is not much larger than the star's (Q*). Both 
of these factors are highly uncertain with various the- 
oretical and observational constraints ranging over sev- 
eral orders of magnitude. In p articular tidal circular- 
ization of main sequence stars dClaret fc Cunhal 119971; 
Meibom fc Mathieul [2005t iZahn fc Bouchetl I1989T IZahnl 



19891) seem to indicate 10 5 < < 10 6 . On the 



other hand, the discovery of extremely short period mas- 



sive planets, the t wo most dramati c being WASP- 18 b 
(IHellier et al.ll2009l) and WASP-19b (IHellier et al.ll20ll . 
seems to be inconsistent with such efficient dissipation 
(Penev et al. in preparation), requiring much larger val- 
ues Q± > 10 s , whic h coincide well wi t h the theoreti- 
cal values derived by iPenev fc Sassel ov (2011), who ar- 
gue that binary stars and star-planet systems are sub- 
ject to different modes of dissipation in the star. The 
tidal dissipation parameter in the planet has also been 
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the subject of ma ny studies attempting to constrain it 
either from theory ( Bodcnh eimer. Laughlin. fc Linll2003t 
lOgilvie fc Linl I2004J ) or from the observed config uration 
of Jupiter's satellites (jGoldreich fc Soteil 119661 ) giving 

10 5 < Qp < 10 7 . 

With this in mind we conclude that the circulariza- 
tion of HAT-P-34b's orbit is likely dominated by the 
tidal dissipation in the planet and using Qp = 10 6 and 
the expression for the tidal c ircularization timescale from 
lAdarns and Laughlin! (|2006[ ). we estimate the eccentricity 
of HAT-P-34b should decay on the scale of 2Gyr, i.e. it is 
not in conflict with theoretical expectations. The possi- 
ble outer companion indicated by the RV trend may also 
be responsible for pum ping the eccentricity of the inner 
planet HAT-P-34b (see lCorreia et al.ll201ll for a discus- 
sion) . 

Figure QT] shows HAT-P-34b on the orbital period- 
eccentricity plane of TEPs with well determined param- 
eters (using our own compilation that attempts to keep 
up with various refinements to these parameters). It is 
apparent that eccentricity is correlated with orbital pe- 
riod and with planet mass, as expected from tidal theory. 
HAT-P-34b lies in a sparse position in these diagrams; 
for example, it has a high eccentricity for its period, the 
only similar planet being HAT-P-2b. 

Fig ure [H is a "tidal" plot (see Fig. 3 of iPont et al.l 
2011), showing TEPs with well measured properties in 
the a/ Rp-Mp/M* plane, using mor e data points (in - 
cluding the present discoveries) than IPont et al.l (|201lD . 

Since t c = (4/63)Q PV /(a 3 jGMl(a/R P ) 5 M p /M ir , we ex- 
pect planets with small relative semi-major axis (a/Rp) 
or planets with small relative mass (M p /M+) to be circu- 
larized. This is indeed the shown by the intensity 
(color) scale representing eccentricity. For hot Jupitcrs 
that migrate in by circularization of an initially very ec- 
centric orbit, the expe cted "parking distance" is ~ 2a# 
(|Ford and Rasidl2006D . where an is the semi-major axis 
at which the radius of the planet equals its Hill radius. 
The thick solid line in Figure [12] shows this relation. A 
fairly good match for the dividing line between the cir- 
cularized (denoted by black points) and eccentric (grey 
or color) points is at a « 4a# (marked with a thin solid 
line). This relation now includes very small mass Kepler 
discoveries. HAT-P-34b belongs to the sparse group of 
high relative semi-major axis (a/R p ) and massive extra- 
solar planets. 
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NNG04GN74G, NNX08AF23G. We acknowledge partial 
funding of the HATNet follow-up effort from NSF AST- 
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A289Hr) and NASA (N167Hr, N029Hr). This paper uses 
observations obtained with facilities of the Las Cumbres 
Observatory Global Telescope. Data presented in this 
paper are based on observations obtained at the HAT 
station at the Submillimeter Array of SAO, and HAT 
station at the Fred Lawrence Whipple Observatory of 
SAO. The authors wish to recognize and acknowledge 
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summit of Mauna Kea has always had within the indige- 
nous Hawaiian community. We are most fortunate to 
have the opportunity to conduct observations from this 
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Fig. 11. — Orbital-period-eccentricity diagram of TEPs with eccentricity uncertainty less than 0.1. The color (greyscale shade) of the 
symbols indicates the mass of each planet. HAT-P-34b is labelled. As expected from tidal evolution theory, high eccentricity planets tend 
to have longer orbital periods and greater masses. 
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Fig. 12. — A "tidal" diagram following Fig. 3 of Pont ct al. (2011). The color (greyscale shade) of the symbols indicates the eccentricity, 
we assume zero eccentricity for planets that have a measured eccentricity within 4cr of zero. The dotted line shows the locus of points with 
a circularization time-scale of 1 Gyr assuming small eccentricity, Qp = 10 6 , and P = 3d. The thick solid line shows the relation a = 2a h 
where ajy is the semi-major axis at which the radius of the planet equals its Hill radius. The thin solid line shows a = 4au- 
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APPENDIX 

SPECTROSCOPIC AND PHOTOMETRIC DATA 

The following tables present the spectroscopic data (radial velocities, bisector spans, and activity index measure- 
ments) and high precision photometric data for the four planets presented in this paper. 



HAT-P-34b, HAT-P-35b, HAT-P-36b, and HAT-P-37b 



TABLE 7 

Relative radial velocities, bisector spans, and activity index measurements 

of HAT-P-34. 



BJD 


RV a 


^RV b 


BS 


CBS 


S c 


Phase 


Instrument 


(2,454,000+) 


(ms- 1 ) 


(ms- 1 ) 


(ms- 1 ) 


(ms -1 ) 








1669.92722 


— 158.62 


50.92 








0.188 


Subaru 


1339.93516 


—201.37 


57.36 








0.190 


Subaru 




—204.06 


44.53 








0.191 


Subaru 


lo41. llooU 


—246.80 


47.11 








0.406 


Subaru 


lo41.12U5o 


on A O *7 

—204.87 


47.32 








0.407 


Subaru 


lo41.124oo 


— 193.24 


43.00 








0.408 


Subaru 


lo i4.12Uoy 


— 196.46 


11.37 


—2.85 


6.57 


0.176 


0.459 


rvcck 


1374.85735 






—9.20 


4.18 


0.173 


0.594 


Keck 


1374.86800 


—80.49 


11.89 


—25.52 


9.34 


0.176 


0.596 


Keck 


1 17IT neon o 

1375.95203 


285.94 


10.24 


— 15.73 


4.46 


0.178 


0.795 


Keck 


1d7o.12o2o 


—265.18 


12.71 


— 13.52 


7.86 


0.174 


0.194 


Keck 


1379.0745s 


—22b. 2b 


12.73 


— 10.85 


3.72 


0.176 


0.368 


Keck 


i o^n ^ f o /"* 

1379.5o5Jo 


— 162.68 


61.60 






' " 


0.458 


FIES 


1 Q Cn 1 Q Q A A 

1doU.1oo44 


— 135.85 


13.51 


60.87 


20.44 


0.180 


0.562 


rveck 


1380.49378 


3.45 


58.70 








0.628 


FIES 


1001 

lool. 11125 


110.87 


14.71 


—9.04 


12.64 


0.175 


0.741 


Keck 


1381.55500 


269.57 


46.50 








0.823 


FIES 


1383.4/949 


— 108.72 


56.00 








0.176 


TTTT7C 

r IrijO 


1384.48968 


-173.76 


50.60 








0.361 


FIES 


1400.85414 


-141.78 


15.43 


27.38 


7.53 


0.167 


0.362 


Keck 


1403.82567 


497.16 


15.18 


38.35 


6.14 


0.173 


0.907 


Keck 


1404.83902 


4.37 


14.05 


0.79 


7.25 


0.174 


0.093 


Keck 


1415.04513 


310.06 


14.03 


2.76 


8.29 


0.168 


0.965 


Keck 


1423.66787 


-115.78 


78.50 








0.546 


FIES 


1424.64318 


41.31 


139.40 








0.725 


FIES 


1425.57489 


474.54 


48.70 








0.896 


FIES 


1426.64513 


-154.35 


57.60 








0.092 


FIES 


1464.85484 


-99.38 


15.03 


-11.41 


11.10 


0.164 


0.100 


Keck 


1465.94874 


-249.34 


13.70 


-19.33 


6.92 


0.163 


0.300 


Keck 


1467.71998 


0.98 


13.61 


-12.70 


8.00 


0.164 


0.625 


Keck 



NOTE. — Note that for the iodine-free template exposures we do not measure the RV but do measure the 
BS and S index. Such template exposures can be distinguished by the missing RV value. 

a The zero-point of these velocities is arbitrary. An overall offset -y re j fitted to these velocities in Scction f^2l 
has not been subtracted. 

Internal errors excluding the component of astrophysical jitter considered in Scction f3"2l 
C Chromosphcric activity index, calibrated to the scale of lVaughan, Preston fc Wilso"nl 1 19781 ). 



TABLE 8 

Relative radial velocities, bisector spans, and activity index measurements 

of HAT-P-35. 



BJD 


RV a 


<mv b 


BS 


OF3S 




s c 


Phase 


Instrument 


(2,454,000+) 


(ms' 1 ) 


(ms" 1 ) 


(ms- 1 ) 


(ms- 1 ) 










1466.12094 


-92.77 


2.40 


-7.06 


3.64 





129 


0.139 


Keck 


1468.11371 


111.39 


1.86 


4.37 


2.92 





129 


0.686 


Keck 


1468.12932 






0.45 


2.76 





129 


0.690 


Keck 


1470.12970 


-118.71 


2.72 


2.26 


2.07 





124 


0.239 


Keck 


1482.73469 d 


106.37 


10.20 


7.70 


9.60 






0.695 


FIES 


1483.74899 d 


39.07 


10.70 


-17.40 


6.50 






0.973 


FIES 


1486.74544 d 


97.37 


8.50 


8.90 


8.30 






0.795 


FIES 


1488.74737 d 


-152.93 


16.90 


24.70 


17.20 






0.344 


FIES 


1490.72460 d 


78.27 


8.50 


-23.80 


6.40 






0.886 


FIES 


1501.03468 


117.01 


2.39 


9.98 


2.63 





127 


0.713 


Keck 


1523.09695 


126.39 


2.09 


0.65 


4.88 





134 


0.763 


Keck 


1529.15460 


-55.78 


2.32 


-11.07 


2.00 





129 


0.424 


Keck 


1545.14866 


103.29 


3.12 


9.13 


3.56 





125 


0.810 


Keck 



NOTE. — Note that for the iodine-free template exposures we do not measure the RV but do measure the 
BS and S index. Such template exposures can be distinguished by the missing RV value. 

a The zero-point of these velocities is arbitrary. An overall offset 7 re j fitted to these velocities in Scction f3"2l 
has not been subtracted. 

k Internal errors excluding the component of astrophysical jitter considered in Scction |3.2| 

C Chromosphcric activity index, calibrated to the scale of |Vaughan, Preston &: WilsoT] j 15781 ). 

d The FIES/NOT observations of HAT-P-35 were not used in the analysis, see the footnote to Table [4] 

Transit ingress began during the hour-long exposure obtained at phase 0.973, and the exposure obtained at 

phase 0.344 has a low S/N ratio and was obtained during morning twilight. 
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TABLE 9 

Relative radial velocities, bisector spans, and activity 
index measurements of hat-p-36. 



BJD 


RV a 


°~RV 


b 


BS 


(JBS 


Phase 


(2,454,000+) 


(m s 


') 


(m 


/ 


fill s 1 


(ms -1 ) 






1547.02139 


-206 


37 


22 


40 


20.50 


15.80 





319 


1550.04236 


200 


83 


22 


80 


-11.40 


21.60 





594 


1554.99394 


-293 


87 


11 


10 


-7.50 


7.50 





325 


1556.05388 


-237 


67 


7 


40 


3.20 


7.70 





123 


1557.05533 


240 


83 


24 


60 


5.10 


18.40 





878 


1558.95741 


-355 


47 


10 


50 


4.30 


8.70 





311 


1559.99613 


-199 


97 


17 


70 


-14.10 


17.40 





093 


1576.04204 


-326 


97 


16 


70 









182 


1576.91790 


331 


53 


14 


10 









842 


1578.04243 


265 


83 


7 


40 









689 


1579.96119 


-277 


57 


9 


00 









135 


1580.92488 


257 


33 


12 


40 









861 



The zero-point of these velocities is arbitrary. An overall offset 7 rcl fitted to 
these velocities in Section |3.2| has not been subtracted. 

Internal errors excluding the component of astrophysical jitter considered in 

Section rr^i 



TABLE 10 

Relative radial velocities, bisector spans, and activity 
index measurements of hat-p-37. 



BJD RV a fT RV b BS cr B s Phase 



(2,454,000+) 


(m s 


X ) 


(m s 


- 1 ) 


(m s 


- 1 ) 


(m s 


- 1 ) 






1646.92799 


146 


43 


35 


50 


26 


20 


16 


30 





710 


1648.01278 


-111 


17 


24 


60 


-28 


00 


16 


30 





098 


1648.91342 


-88 


77 


23 


20 





60 


22 


20 





420 


1649.89870 


172 


53 


29 


10 


8 


30 


24 


10 





772 


1650.96241 


-102 


97 


17 


50 


12 


30 


16 


90 





153 


1652.92638 


128 


33 


26 


30 


-13 


70 


15 


00 





855 


1656.98738 


-238 


57 


12 


80 


-28 


70 


23 


90 





306 


1658.96338 


-19 


87 


20 


90 


-9 


00 


17 


00 





013 


1663.90569 


169 


23 


20 


60 


-5 


60 


15 


10 





779 


1664.91140 


-163 


77 


25 


10 


-28 


30 


19 


50 





139 


1665.91961 


62 


93 


32 


60 


-29 


50 


16 


80 





499 


1671.84933 


125 


23 


35 


30 


47 


40 


27 


20 





619 


1698.79464 


-143 


97 


25 


00 


48 


00 


16 


20 





251 



The zero-point of these velocities is arbitrary. An overall offset 7 rcl fitted to 
these velocities in Scction fT^ has not been subtracted. 

Internal errors excluding the component of astrophysical jitter considered in 
Section [3"2l 



TABLE 11 

High-precision differential photometry of HAT-P-34. 



BJD Mag a cr M ag Mag(orig) b Filter 



(2,400,000+) 



55338.79359 


0.00011 


55338.79427 


-0.00108 


55338.79494 


0.00263 


55338.79563 


-0.00041 


55338.79633 


0.00158 


55338.79701 


0.00036 


55338.79769 


-0.00243 


55338.79837 


-0.00305 


55338.79905 


0.00184 


55338.79969 


0.00187 



0.00080 10.12720 i 

0.00080 10.12640 i 

0.00080 10.12910 i 

0.00080 10.12660 i 

0.00080 10.12850 i 

0.00079 10.12720 i 

0.00079 10.12440 i 

0.00079 10.12340 i 

0.00079 10.12940 i 

0.00079 10.12900 i 



NOTE. — This table is available in a machine-readable form in the online journal. 
A portion is shown here for guidance regarding its form and content, 

a The out-of-transit level has been subtracted. These magnitudes have been 
subjected to the EPD and TFA procedures, carried out simultaneously with the 
transit fit. 

Raw magnitude values without application of the EPD and TFA procedures. 



HAT-P-34b, HAT-P-35b, HAT-P-36b, and HAT-P-37b 



TABLE 12 

High-precision differential photometry of HAT-P-35. 



BJD 

(2,400,000+) 


Mag a 


"Mag 


Magi 


orig) b 


55578.68612 


0.00996 





00115 


11 


35390 


55578.68766 


0.01222 





00104 


11 


35670 


55578.68921 


0.01205 





00101 


11 


35720 


55578.69075 


0.00983 





00097 


11 


35550 


55578.69228 


0.01193 





00099 


11 


35670 


55578.69383 


0.01274 





00113 


11 


35760 


55578.69504 


0.01155 





00113 


11 


35620 


55578.69622 


0.01243 





00115 


11 


35790 


55578.69742 


0.01049 





00114 


11 


35540 


55578.69863 


0.00960 





00114 


11 


35430 



Filter 



NOTE. — This tabic is available in a machine- readable form in the online 
journal. A portion is shown here for guidance regarding its form and content. 

a The out-of-transit level has been subtracted. These magnitudes have been 
subjected to the EPD and TFA procedures, carried out simultaneously with the 
transit fit. 

k Raw magnitude values without application of the EPD and TFA procedures. 



TABLE 13 

High-precision differential photometry of HAT-P-36. 



BJD 

(2,400,000+) 



Mag a 



"Mag Mag(orig) b Filter 



55555 


84870 





00292 





00081 


10 


95490 


55555 


85050 





00957 





00083 


10 


96050 


55555 


85205 





00817 





00078 


10 


96100 


55555 


85362 





01416 





00079 


10 


96370 


55555 


85516 





01203 





00079 


10 


96260 


55555 


85671 





01521 





00079 


10 


96600 


55555 


85826 





01690 





00078 


10 


96760 


55555 


85981 





01833 





00079 


10 


96870 


55555 


86135 





01590 





00078 


10 


96410 


55555 


86292 





01733 





00078 


10 


96540 



NOTE, — This table is available in a machine-readable form in the online 
journal. A portion is shown here for guidance regarding its form and content. 

a The out-of-transit level has been subtracted. These magnitudes have been 
subjected to the EPD and TFA procedures, carried out simultaneously with the 
transit fit. 

k Raw magnitude values without application of the EPD and TFA procedures. 



TABLE 14 

High-precision differential photometry of HAT-P-37. 



BJD 


Mag a 


"Mag 


Mag(orig) b Filter 


(2,400,000+) 








55616.96855 


0.02527 


0.00133 


12.41350 i 


55616.97081 


0.02387 


0.00135 


12.41200 i 


55616.97308 


0.02220 


0.00133 


12.41020 i 


55616.97540 


0.02311 


0.00131 


12.41290 i 


55616.97774 


0.02297 


0.00132 


12.41110 i 


55616.97986 


0.02149 


0.00133 


12.40940 i 


55616.98176 


0.02186 


0.00133 


12.40980 i 


55616.98366 


0.01909 


0.00131 


12.40740 i 


55616.98555 


0.02301 


0.00132 


12.41180 i 


55616.98744 


0.02165 


0.00129 


12.40910 i 



NOTE. — This table is available in a machine-readable form in the online 
journal. A portion is shown here for guidance regarding its form and content. 

a The out-of-transit level has been subtracted. These magnitudes have been 
subjected to the EPD and TFA procedures, carried out simultaneously with the 
transit fit. 

k Raw magnitude values without application of the EPD and TFA procedures. 



